Development of an Advanced Force Field for Water using Variational Energy Decomposition Analysis

Given the piecewise approach to modeling intermolecular interactions for force fields, they can be difficult to parameterize since they are fit to data like total energies that only indirectly connect to their separable functional forms. Furthermore, by neglecting certain types of molecular interactions such as charge penetration and charge transfer, most classical force fields must rely on, but do not always demonstrate, how cancellation of errors occurs among the remaining molecular interactions accounted for such as exchange repulsion, electrostatics, and polarization. In this work we present the first generation of the (many-body) MB-UCB force field that explicitly accounts for the decomposed molecular interactions commensurate with a variational energy decomposition analysis, including charge transfer, with force field design choices that reduce the computational expense of the MB-UCB potential while remaining accurate. We optimize parameters using only single water molecule and water cluster data up through pentamers, with no fitting to condensed phase data, and we demonstrate that high accuracy is maintained when the force field is subsequently validated against conformational energies of larger water cluster data sets, radial distribution functions of the liquid phase, and the temperature dependence of thermodynamic and transport water properties. We conclude that MB-UCB is comparable in performance to MB-Pol, but is less expensive and more transferable by eliminating the need to represent short-ranged interactions through large parameter fits to high order polynomials

The Poisson-Boltzmann model for implicit solvation of electrolyte solutions: Quantum chemical implementation and assessment via Sechenov coefficients.

We present the theory and implementation of a Poisson-Boltzmann implicit solvation model for electrolyte solutions. This model can be combined with arbitrary electronic structure methods that provide an accurate charge density of the solute. A hierarchy of approximations for this model includes a linear approximation for weak electrostatic potentials, finite size of the mobile electrolyte ions, and a Stern-layer correction. Recasting the Poisson-Boltzmann equations into Euler-Lagrange equations then significantly simplifies the derivation of the free energy of solvation for these approximate models. The parameters of the model are either fit directly to experimental observables-e.g., the finite ion size-or optimized for agreement with experimental results. Experimental data for this optimization are available in the form of Sechenov coefficients that describe the linear dependence of the salting-out effect of solutes with respect to the electrolyte concentration. In the final part, we rationalize the qualitative disagreement of the finite ion size modification to the Poisson-Boltzmann model with experimental observations by taking into account the electrolyte concentration dependence of the Stern layer. A route toward a revised model that captures the experimental observations while including the finite ion size effects is then outlined. This implementation paves the way for the study of electrochemical and electrocatalytic processes of molecules and cluster models with accurate electronic structure methods

Trap models with slowly decorrelating observables

We study the correlation and response dynamics of trap models of glassy dynamics, considering observables that only partially decorrelate with every jump. This is inspired by recent work on a microscopic realization of such models, which found strikingly simple linear out-of-equilibrium fluctuation-dissipation relations in the limit of slow decorrelation. For the Barrat-Mezard model with its entropic barriers we obtain exact results at zero temperature $T$ for arbitrary decorrelation factor $\kappa$. These are then extended to nonzero $T$, where the qualitative scaling behaviour and all scaling exponents can still be found analytically. Unexpectedly, the choice of transition rates (Glauber versus Metropolis) affects not just prefactors but also some exponents. In the limit of slow decorrelation even complete scaling functions are accessible in closed form. The results show that slowly decorrelating observables detect persistently slow out-of-equilibrium dynamics, as opposed to intermittent behaviour punctuated by excursions into fast, effectively equilibrated states.Comment: 29 pages, IOP styl

Cardiovascular disease in a cohort exposed to the 1940-45 Channel Islands occupation

BACKGROUND To clarify the nature of the relationship between food deprivation/undernutrition during pre- and postnatal development and cardiovascular disease (CVD) in later life, this study examined the relationship between birth weight (as a marker of prenatal nutrition) and the incidence of hospital admissions for CVD from 1997–2005 amongst 873 Guernsey islanders (born in 1923–1937), 225 of whom had been exposed to food deprivation as children, adolescents or young adults (i.e. postnatal undernutrition) during the 1940–45 German occupation of the Channel Islands, and 648 of whom had left or been evacuated from the islands before the occupation began. METHODS Three sets of Cox regression models were used to investigate (A) the relationship between birth weight and CVD, (B) the relationship between postnatal exposure to the occupation and CVD and (C) any interaction between birth weight, postnatal exposure to the occupation and CVD. These models also tested for any interactions between birth weight and sex, and postnatal exposure to the occupation and parish of residence at birth (as a marker of parish residence during the occupation and related variation in the severity of food deprivation). RESULTS The first set of models (A) found no relationship between birth weight and CVD even after adjustment for potential confounders (hazard ratio (HR) per kg increase in birth weight: 1.12; 95% confidence intervals (CI): 0.70 – 1.78), and there was no significant interaction between birth weight and sex (p = 0.60). The second set of models (B) found a significant relationship between postnatal exposure to the occupation and CVD after adjustment for potential confounders (HR for exposed vs. unexposed group: 2.52; 95% CI: 1.54 – 4.13), as well as a significant interaction between postnatal exposure to the occupation and parish of residence at birth (p = 0.01), such that those born in urban parishes (where food deprivation was worst) had a greater HR for CVD than those born in rural parishes. The third model (C) found no interaction between birth weight and exposure to the occupation (p = 0.43). CONCLUSION These findings suggest that the levels of postnatal undernutrition experienced by children, adolescents and young adults exposed to food deprivation during the 1940–45 occupation of the Channel Islands were a more important determinant of CVD in later life than the levels of prenatal undernutrition experienced in utero prior to the occupatio

Fluvial Volumes, Timescales, and Intermittency in Milna Crater, Mars

Ancient lake deposits and valley networks on Mars provide strong evidence that its surface was once modified by liquid water, but the extent of that modification is still debated. Ancient lacustrine deposits in Milna Crater provide insight into the timescale and fluid volume required to construct fluvially derived sedimentary deposits near the Noachian-Hesperian boundary. Placing the lacustrine deposits their regional context in Paran Valles provides a quantitative measurement of the intermittency of large, water-mediated sediment transport events in that region

Response of Methanogens in Arctic Sediments to Temperature and Methanogenic Substrate Availability

Although cold environments are major contributors to global biogeochemical cycles, comparatively little is known about their microbial community function, structure, and limits of activity. In this study a microcosm based approach was used to investigate the effects of temperature, and methanogenic substrate amendment, (acetate, methanol and H2/CO2) on methanogen activity and methanogen community structure in high Arctic wetlands (Solvatnet and Stuphallet, Svalbard). Methane production was not detected in Stuphallet sediment microcosms (over a 150 day period) and occurred within Solvatnet sediments microcosms (within 24 hours) at temperatures from 5 to 40°C, the maximum temperature being at far higher than in situ maximum temperatures (which range from air temperatures of -1.4 to 14.1°C during summer months). Distinct responses were observed in the Solvatnet methanogen community under different short term incubation conditions. Specifically, different communities were selected at higher and lower temperatures. At lower temperatures (5°C) addition of exogenous substrates (acetate, methanol or H2/CO2) had no stimulatory effect on the rate of methanogenesis or on methanogen community structure. The community in these incubations was dominated by members of the Methanoregulaceae/WCHA2-08 family-level group, which were most similar to the psychrotolerant hydrogenotrophic methanogen Methanosphaerula palustris strain E1-9c. In contrast, at higher temperatures, substrate amendment enhanced methane production in H2/CO2 amended microcosms, and played a clear role in structuring methanogen communities. Specifically, at 30°C members of the Methanoregulaceae/WCHA2-08 predominated following incubation with H2/CO2, and Methanosarcinaceaeand Methanosaetaceae were enriched in response to acetate addition. These results may indicate that in transiently cold environments, methanogen communities can rapidly respond to moderate short term increases in temperature, but not necessarily to the seasonal release of previously frozen organic carbon from thawing permafrost soils. However, as temperatures increase such inputs of carbon will likely have a greater influence on methane production and methanogen community structure. Understanding the action and limitations of anaerobic microorganisms within cold environments may provide information which can be used in defining region-specific differences in the microbial processes; which ultimately control methane flux to the atmosphere

Distribution of thermophilic endospores in a temperate estuary indicate that dispersal history structures sediment microbial communities

Endospores of thermophilic bacteria are found in cold and temperate sediments where they persist in a dormant state. As inactive endospores that cannot grow at the low ambient temperatures, they are akin to tracer particles in cold sediments, unaffected by factors normally governing microbial biogeography (e.g., selection, drift, mutation). This makes thermophilic endospores ideal model organisms for studying microbial biogeography since their spatial distribution can be directly related to their dispersal history. To assess dispersal histories of estuarine bacteria, thermophilic endospores were enriched from sediments along a freshwater‐to‐marine transect of the River Tyne in high temperature incubations (50°C). Dispersal histories for 75 different taxa indicated that the majority of estuarine endospores were of terrestrial origin; most closely related to bacteria from warm habitats associated with industrial activity. A subset of the taxa detected were marine derived, with close relatives from hot deep marine biosphere habitats. These patterns are consistent with the sources of sediment in the River Tyne being predominantly terrestrial in origin. The results point to microbial communities in estuarine and marine sediments being structured by bi‐directional currents, terrestrial run‐off and industrial effluent as vectors of passive dispersal and immigration

Socioeconomic Inequalities in Disability-free Life Expectancy in Older People from England and the United States: A Cross-national Population-Based Study

BACKGROUND: We examined socioeconomic inequalities in disability-free life expectancy in older men and women from England and the United States and explored whether people in England can expect to live longer and healthier lives than those in the United States. METHODS: We used harmonized data from the Gateway to Global Aging Data on 14,803 individuals aged 50+ from the U.S. Health and Retirement Study (HRS) and 10,754 from the English Longitudinal Study of Ageing (ELSA). Disability was measured in terms of impaired activities and instrumental activities of daily living. We used discrete-time multistate life table models to estimate total life expectancy and life expectancy free of disability. RESULTS: Socioeconomic inequalities in disability-free life expectancy were of a similar magnitude (in absolute terms) in England and the United States. The socioeconomic disadvantage in disability-free life expectancy was largest for wealth, in both countries: people in the poorest group could expect to live seven to nine fewer years without disability than those in the richest group at the age of 50. CONCLUSIONS: Inequalities in healthy life expectancy exist in both countries and are of similar magnitude. In both countries, efforts in reducing health inequalities should target people from disadvantaged socioeconomic groups

Nonlinear porous medium flow with fractional potential pressure

We study a porous medium equation, with nonlocal diffusion effects given by an inverse fractional Laplacian operator. We pose the problem in n-dimensional space for all t>0 with bounded and compactly supported initial data, and prove existence of a weak and bounded solution that propagates with finite speed, a property that is nor shared by other fractional diffusion models.Comment: 32 pages, Late

Effect of spatial bias on the nonequilibrium phase transition in a system of coagulating and fragmenting particles

We examine the effect of spatial bias on a nonequilibrium system in which masses on a lattice evolve through the elementary moves of diffusion, coagulation and fragmentation. When there is no preferred directionality in the motion of the masses, the model is known to exhibit a nonequilibrium phase transition between two different types of steady states, in all dimensions. We show analytically that introducing a preferred direction in the motion of the masses inhibits the occurrence of the phase transition in one dimension, in the thermodynamic limit. A finite size system, however, continues to show a signature of the original transition, and we characterize the finite size scaling implications of this. Our analysis is supported by numerical simulations. In two dimensions, bias is shown to be irrelevant.Comment: 7 pages, 7 figures, revte